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^ ABSTRACT 

O 
O 

CN ! The 3+1 Hamiltonian formulation in the gauge DfN = —K on the lapse 

'i> ' function fixes the direction of time associated with the trace K of the extrinsic 



O 



curvature tensor. The Hamiltonian equations hereby become hyperbolic. We 
study this new system for black hole spacetimes that are asymptotically quies- 
cent, which introduces analyticity properties that can be exploited for numerical 
O ' I calculations by compactification in spherical coordinates with complex radius fol- 

lowing a Mobius transformation. Conformal flat initial data of two black holes 



(3JQ| are hereby invariant, and correspond to a turn point in a pendulum, up for a pair 
of separated black holes and down for a single black hole. Here, Newton's law 

^ _ appears in the relaxation of / = 2 deformations of semi-infinite poloidal surface 

^ ' elements, defined by the moment of inertia of the binary. 
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Introduction 



The calculation of gravitational radiation produced in the merger of two black holes is of 
considerable theoretical interest. It may also provide templates of wave- forms for analysis of 
^ ' LIGO- Virgo data. To be of practical interest to matched filtering, the calculational methods 
I will have to recover many wave-periods over an extended parameter range, such as in the 
problem of binary black hole coalescence with a range of black hole masses and spins. It 
poses the challenge of efficient and stable phase-accurate numerical methods. It has been 
appreciated that this requires an inherently stable formulation of the dynamical evolution 
of spacetime. Spectral methods provide the most efficient representation of functions that 
are periodic and everywhere analytic, and preserve accurate phase information in hyperbolic 
evolution. 

Here, we consider the 3+1 Hamilt onian equations for the dynamical behavior of space- 



time (lArnowitt. Deser fc Misneiill962l ). They are attractive, in describing the evolution in 
terms of quantities that have unambiguous geometric interpretations. However, for arbi- 
trary choices of gauge, they are known to give rise to computational instabilities, which have 
thusfar prevented their application to large-scale numerical relativity. 



- 2 - 



Here, we propose a new gauge for time-evolution, which describes a correlation between 
the lapse function and the extrinsic curvature of the foliation of three-surfaces in space- 
time. Linearized stability analysis shows that the full system of Hamiltonian equations now 
becomes hyperbolic with respect to arbitrary perturbations of the three-metric. Ne xt, we 



focus on black hole spacetimes which are asymptotically quiescent (Ivan PuttenI 120061 ). This 
asymptotic property, if present in the initial data, is preserved for all time by causality in 
the equations of general relativity. It can be exploited by compactification in spherical coor- 
dinates with complex radius, {z, 9, 0) with complex radius z, — oo + is < z < oo + is, where 
E/s denotes a dimensional continuation parameter for a spacetime with total mass energy E. 
Here, s can be chosen sufficiently large, i.e., at least on the order of E, to avoid singularities 
associated with black holes, even those that may carry spin. In this approach, we focus on 
the outer expansion of the Green's function. It makes possible a three-dimensional spectral 
representations after application of a Mobius transformation (which preserves conformally 
flat initial data). 

The combination of hyperbolic Hamiltonian equations and spectral representations offers 
a novel starting point for studying the dynamical evolution of black hole spacetimes. For 
illustrative purposes, we consider the conformally flat initial data of Schwarzschild black 
holes, as in the problem of a head-on collision of two black holes. Here, conformally flat 
data are special, in representing turnings points corresponding to the up or down position 
of a pendulum. The proposed lapse function is natural choice about such turning points. To 
illustrate compactiflcation by a Mobius transformation, we explicitly identify Newton's law 
in terms of the evolution of the quadrupole deformations in area of poloidal flux surfaces, 
given by the moment of inertia / of spacetime. 



2. Asymptotic wave motion 

At large distances from the source region, the two wave-modes in the gravitational 
fleld appear as planar waves in the three-metric hij, commonly expressed in the transverse 
traceless gauge as 

/I \ 

hij=\0 1 + b c , b = Bsm{kz), c = Csm{kz) (1) 
\ c 1-b J 

B, C denote the amplitude of -|- and x plane-waves along the direction with wave-number 
k. The three-dimensional Ricci tensor associated with hij hereby assumes the following 
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structure. For the + waves with {B,C) = (5,0), we have 

JDK * ^ 

R+ = ^ [ - sin A;^ - cos^ kz + sin^ kz | , (2) 

sin kz — B cos^ kz — B^ sin^ kz 



B^ 
'2A 



where b = B sin kz, c = and A = 1 — 5^ sin^ kz, and for the x waves with {B, C) = (0, C), 
b = 0, c = C sin kz, we find, similarly, 

= ^ cos2 kz -A sin A;^ | , (3) 
\ —A sin kz —C cos^ kz 

where A = 1 — sin^ kz. Hence, we have the following exact algebraic identities: 

Rt2-Rtz = -Bk^Mkz), i?+3 = 0, (4) 
^3'3-^22 = 0, R^^^-^esin{kz), (5) 

where the numerical indices refer to normalized directions in the coordinate directions 
(r,^,0). 

The identies ([5]) show that the two polarization modes satisfy traveling wave solutions 

b{t,r) = B smrj, c{t,r) = C smrj, t] = k(t — r) (6) 

in the large-distance approximation (neglecting 1/r curvature terms, in this section) with 
vanishing lapse functions /3i = and constant lapse function = 1, satisfying 

{h22 - h,s) = -2dt{K,, - K22) = 2{R22 - R33), 5>23 = -29^7^23 = 2i?23, (7) 

where b{t,r) = | (/122 — ^11) , c(t,r) = hi2, using the numerical indices to denote the 
normalized directions in (r, 6, 0). The angular directions "22-33" and "23" therefore represent 
the hyperbolic directions, which contain all wave-motion. 



3. General 3-|-l decomposition of the metric 



W e consider the 3+1 decomposition of the line-element in terms of hij (IThorne. Price &: McDonald 
1986} ) 



ds"^ = -N^dt^ + hij {dx' + (3'dt){dx^ + (3jdt) 
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The lapse and shift functions can be seen to define the motion of observers with velocity four- 



vector = (1,0,0,0), for example, by considering their acceleration a* 



dt 



U 



-r 



00 



~\g'' {2gio - goo,i) = -h'' (^A - NdiN + (3Wi(3jj = ^-^NdiN, which is determined entirely 

by the redshift factor when (3i = j3i = 0. In general, the lapse and shift functions appear 
in the conservation of energy and momentum, of test particles moving along geodesies. We 
need not impose the time-symmetric gauge condition with /5j = 0. The three-metric hij 
satisfies the Hamiltonian evolution equations 



(9) 



where DtKij = 
three-tensor W^, 



dtKij 



Rij + N-^DiDjN. 



KjmDjP"^ and Wij denotes the gauged 



We propose the time-evolution accoring to 

= -K, 




-2KijN, 



(10) 



where 



Rij and /C = KKij 



2KY'K 



The gauge condition DtN = —K is curvature- driven evolution. It is different from the 
product of curva ture and lapse fu n ction in the ha rmonic s licing condition dtN = —N^K in 
Eqs. (69)-(77) of [Abrahams et al.l ( 19971 ): see also [Brown J2008h for a recent review. 



In the asymptotically flat region about hi 



6ij and = 1, we have 



dtN = -K, d^hij = -2Rij + 2DiDjN, d^K = AK. 



We recall that (IWaldlll984h 



1 1 - 1 - 

--A6hij + +-did''6hej + -djd^Shei 



(12) 



where hij = 6hij — ^6ij6g, where 6h = h^^dhij refers to the trace of the metric perturbations. 



For harmonic perturbations of the for 5hij ~ hijC *'^*e*'^'^\ we have, with conservation of 



momentum, D^K. 



-iujN 



-K, 6hi 



-2iuj ^Kij, k'Kij 



kjK, so that 



did^h 



lU) 



-1/ 



'2kik^Kej + kikjK) = —ico kikjK 



(13) 



and hence 



R, 



V 



didjN = -k'^hij — iu ^kikjK + iuo ^kikjK = -k^hij 



(14) 
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We conclude that dthij = —2Rij + 2DiDjN gives rise to the dispersion relation 

00^ = e (15) 

for arbitrary, small amplitude metric perturbations. We conclude that the system K in flTUl) 
is asymptotically stable by virtue of the gauge choice DfN = —K. This contrasts with the 
more common derivation of f|T5|) in the so-called transverse traceless gauge, or harmonic 
coordinates-neither of these two coordinate conditions are used here. 



Conformal fiat initial data 



The conformal decomposition of the Ricci tensor for hij = (p'^^Qij satisfies ( lWaldlll984l ) 

^''^R.jih) R^j{g) - m(f)-^[{n - 2)DiDj4> + gijAcf)] + m^-^ [(1 + m){n- 2)Di(f)Dj(f) + (1 - m * 
{n - 2))gijDP(j)Dp(j)]. For n = 3 and m = 2, it reduces to ^^^Rij =^^^ Rij{g) - 2(j)-^[DiDj(j) + 
gijA(j)] + 2(j)~'^[3Di(j)Dj(j) — gijD'^ipDpCf)], giving the familiar result 

Rh = [Rg - 80-iA,0] . (16) 
Time-symmetric initial data around one or multiple non-rotating black holes can be 



conveniently described in a conformally flat representation (jLindquistlll963l : lBrill fc Lindquist 



19641 : 1 Jansen et al.ll2003l ) on the basis of f |T6l) . It gives rise to solutions for vacuum spacetimes 
in terms of the Green's function G(x*,p*) = Gp of flat spacetime, where denotes the 
position of a point source. 

Spacetime of a single black hole is given by the conformal factor 

$ = 1 + 2TxmGp = 1 + § + ^^i(^) + ^^2(x) + • • • (17) 

where E = m, I = mp^ and the expansion refers to the outer expansion of Gq in complex 
spherical coordinates (z, x, 0), x = cos 9, by exploiting analyticity at infinity, i.e., asymptotic 
quiescence. The same construction applies to two black holes, 

$ = 1 + 27iMGg + 27imGp = 1 + ^ + ^^^^^P,{x) + ^P2{x) + ■■■ (18) 

where E = M + m and / = Mq"^ + mp"^. Without loss of generality, we may choose 
Mq + mp = 0. According to the Hamiltonian evolution equations, the associated shift 
function, for a static foliation of spacetime, i.e., DtKij = 0, satisfies 

AhN + 3NR = 0. (19) 



-6- 




Fig. 1. — The interaction of two Schwarzschild black holes is attractive as seen in the space of 
the complex radial coordinate z (solid curved line) or equivalently repulsive, following closure over 
infinity in the domain of convergence of the outer expansion of the Green's function, i.e., in the 
space of the transformed variable w = Xjz (dotted curved lines). For the latter, Newton's law 
represents the leading order evolution of the area of the poloidal surface elements \z\ > p, governed 
by the moment of inertia associated with the separation p between the two black holes. 
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For the conformally flat data at hand, it reduces to (^'^diN)i = 0, or A?" = (2 - $)/$. 

The leading-order structure of spacetime can be studied by considering poloidal surface 
elements dA — }iJ^}iJzdQdz in the outer region \z\ > p. A measure for the I — 2 deformations 
in surface elements is 

A2{z,t)^ / dA^ / ^^zdedz + --- (20) 

Jz Jo J z Jo 

in z > where ■ ■ ■ refers to time-independent divergences. It gives a geometric equivalence 
about w = {w = 1/z) to the moment of intertia. 



5. Newton's law in the initial quadrupole evolution 

The gauged tensor Wij = —Rij + N^^DiDjN has a regular Taylor series expansion 
about infinity, for the given asymptotically quiescent initial data. For the equal mass, two- 
hole solution with total mass-energy E = 2M, moment of intertia / = 2Mp^, and symmetric 
position about the origin, the conformal factor 



we find the large- 2; asymptotics 
EI 



W* -- 

2z^ 



I 2P2 -3zPi 

-SzPi -z-'[2P2-Po]/{l-x^) |, (22) 

\ -z'^Po{l-x^) 



satisfying W* = with 

W*:W*^ E^P 0^-12 - 9Ez-'^ + (^^E^ - QE'^Ix^^ z'^"^^ + 0{z-'') (23) 

The time-evolution of the / = 2 deformation of poloidal surface elements is given by 

/*7r POO 

A2{z,t)= / dl<lei{hijf/^dedz (24) 

^0 J z 

where I, J — z, 9, where det(/i/j) = heohzz — {hzoY- For the time-symmetric initial data, the 
quadratic off-diagonal terms make no initial contribution, leaving 

A2{z,t)^ / dUhli'hli')dedz^- / {h-i/\li%e + h-,'/'hli%,)d9dz. (25) 

^0 Jz ^ Jo Jz 
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At t = 0, the integrand reduces in the gauge K = with the additional time- symmetric 
gauge condition /3j = at t = (a Newtonian gauge choice, as in a turning point) to 

^ (/i,, + z~%g) z = = ^-2 OK^^ = -§^ + 0{z~'). (26) 

For completely time-symmetric initial data, we thus recover Newton's law for two equal mass 
black holes of mass M and separation 2p in A2{z,t) = —}^EIz~^ + 0{z^^), we find 

M , , 

for initial evolution of A2 at the maximal extent z = p (the lower limit of z) of the poloidal 
surface elements (in the outer expansion considered). 

Thus, for a fixed 1^1 > p, it takes t = 0{z^) to relax Iz^^ by the leading-order z~^ forcing 
terms in ( l22l) . when applied to the asymptotic equations. The asymptotics Iz~^ for large z 
hereby persists for all finite time. This observation is consistent with the a priori notion of 
preserving asymptotic quiescence. Here, Newton's law (p71) associated with the deformations 
of the poloidal surface elements 1^1 > p represent the underlying non-local interaction. 



6. Conclusions 

We have shown that the 3+1 Hamiltonian equations become hyperbolic, when the evolu- 
tion of the lapse function is driven by the extrinsic curvature tensor. The resulting hyperbolic 
system forms a starting point for three-dimensional spectral methods in numerical relativity, 
after compactification by a Mobius transformation for spacetimes that satisfy asymptotic 
quiescence, wherein Newton's law arises in the dynamics of the quadrupole deformations of 
poloidal surface elements. Future work will focus on some numerical examples. 
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